function y = computeTractionLaw(element,etilde) %compute g
ef = element.ef;
ft = element.ft;

%exponential law
y = ft*exp(-etilde/ef);

%linear law
% if etilde < ef
%     y = ttilde*(1 - etilde/ef);
% else
%     y = 0;
% end

end
